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We study the SU(2)k Wess-Zumino-Novikov-Witten (WZNW) theory perturbed by the 
trace of the primary field in the adjoint representation, a theory governing the low-energy be¬ 
haviour of a class of strongly correlated electronic systems. While the model is non-integrable, 
its dynamics can be investigated using the numerical technique of the truncated confor¬ 
mal spectrum approach combined with numerical and analytical renormalization groups 
(TCSA+RG). The numerical results so obtained provide support for a semiclassical analysis 
valid at k 1. Namely, we find that the low energy behavior is sensitive to the sign of the 
coupling constant, A. Moreover for A > 0 this behavior depends on whether k is even or 
odd. With k even, we find definitive evidence that the model at low energies is equivalent 
to the massive 0(3) sigma model. For k odd, the numerical evidence is more equivocal, but 
we find indications that the low energy effective theory is critical. 


I. INTRODUCTION 


Conformal field theories (CFT) describe universal critical behavior and by virtue of this play an 
enormously important role in the physics of strongly correlated systems. This universality is not 
completely lost in the presence of perturbations since, as a rule, the number of relevant operators 
is finite and once restricted by symmetry, often number but a few. Physics of perturbed critical 
models can be rich and complex, especially when the perturbation is non-integrable. For examples 


one may look at the quantum Ising model perturbed simultaneous 
field and the thermal operator [l, Q] or double sine-Gordon models 


by a longitudinal magnetic 


iii 


A focus on relevant perturbations of a CFT is most appropriate when the perturbations are 
strongly relevant. Indeed, the more relevant the perturbation the smaller is the energy scale over 
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which the spectrum is significantly altered. This feature lies at the foundation of the truncated 


conformal spectrum approach (TCSA) introduced in [5]. In the simplest version of this approach 
(TCSA), one truncates the spectrum of the unperturbed CFT which reduces the problem to nu¬ 
merical diagonalization of finite size matrices. Later this idea was combined with a numerical 


renormalization group 


(TCSA+NRG). The TCSA+NRG has been used to tackle a number 

a la], to 


of problems ranging from the excitonic spectrum in semiconducting carbon nanotubes 
quenches in the Lieb-Liniger model [j], kJ , to studying theories whose fields live on a non-compact 


manifold jll]. In a further development, the precision of TCSA or TCSA+NRG computations 


can be improved further upon using perturbative renormalization group techniques 


B BE 


-151]. 


These same renormalization group techniques allow one to use the TCSA to predict gaps in actual 
material systems which possess a finite bandwidth/cutoff [§]. 

Below we will apply the TCSA+NRG to study the (l+l)-dimensional SU(2)k Wess-Zumino- 
Novikov-Witten (WZNW) model perturbed by the trace of the adjoint operator. This is a strongly 
relevant operator with scaling dimension d = ideal for application of the TCSA+NRG method. 


Q 


This perturbed conformal field theory appears in applications such as theories of spin ladders 
(see also the Appendix). A variant of this theory, perturbing SU( 2)*, by the trace of the adjoint on 
the boundary of the system, describes a particular class of Kondo models [17]. Another variant of 
the model, with an additional current-current perturbation, a ppe ared in the description of fermionic 
cold atoms loaded into a one-dimensional optical lattice [18], il9| • 

The perturbed CFT is not integrable except at k = 2 when it is equivalent to the theory of 
three massive Major ana fermions. Below we will use the semiclassical approximation to analyze 
the case of k 1 while using TCSA+NRG for small finite values of k. Our investigations yield 
concrete predictions for the vacuum structure and low-energy excitations for systems described by 
this perturbed conformal field theory. 

The most striking property of the theory is the dependence of its properties on the sign of the 
coupling constant, A. For A > 0 there is a dichotomy in behavior between even and odd k. For odd 
k the semiclassical analysis predicts a massless RG flow from the SU(2)k critical point to SU(2)\ 
in the infrared. The spectrum for even k is always massive and the lowest multiplet is a triplet. 
However the size of the mass depends on the sign of A so that the mass is smaller for A > 0 with 
the ratio m{— A)/m(A) increasing exponentially with k. 
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II. THE PERTURBED SU(2) k WZNW MODEL 


The model in which we are interested is the SU( 2)*, WZNW model perturbed by the trace of 
the primary field in the adjoint representation: 

3 

S = kW(g) + A]T 

a —1 

W(g) = J d 2 xTr(d^g ] d fl g) + T(g); 

T (s) = 2iyr j R ^'Tr^d a g^dpgg^d^g), (2) 

where r(g) is the famous Wess-Zumino term and cr a are the Pauli matrices. The perturbation is 
equivalent to the trace of the WZNW principal field in the adjoint representation: 

TVKVV] ~ (3) 

This is a strongly relevant operator with scaling dimension d = and as such it generates a 
characteristic energy scale 

m ~ | A| 1 /( 2 — 5 (4) 


d 2 xTi-[a a ga a g H 


(1) 


below which the spectrum is strongly modified. 

It is interesting to note that the model perturbed by a single component of matrix & a dj is 


integrable. It was demonstrated in 


2l| that the perturbed Hamiltonian decomposes into a massless 


U(l) CFT and a massive Zj, CFT per turbed by the thermal operator. The properties of the latter 
massive theory were studied in [21,|22j. However with the inclusion of the entire trace, integrability 
is lost for k > 2. At k = 2 the model is equivalent to the model of three massive Major ana fermions 
with mass m ~ A. In this form it has been used to describe the spin 5 = 1 spin ladder fl¬ 
it is interesting to note that for k = 4 when the central charge of the critical WZNW theory is 
equal to c = 2, the model can be recast in abelian form: 

3 


/' 


S = d 2 x 


1 

-87T 


= 


= 2/3, 


( 5 /A) 2 - a y 

a= 1,2 i= 1 

eWeO’A = - 


cos e. 


( 5 ) 


= -1/3, 


as well as the SU( 3)i WZNW model perturbed by the trace of the matrix operator: 


S = W[SU(3y,g] + \Tr(g + g'). 


(6) 
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To get a qualitative understanding of the spectrum of model (Eqn. [TJ), we consider the case 
1 where the model can be treated semiclassically. Using the identity 


^Tr[a a ga a g ] ] = 2 TrgTrg ] - 2 

a= 1 

and the fact that the SU(2) matrix g can be written as 

g = n$i + icr a n a , rig + n 2 = 1, 
we obtain the perturbation in the form 


E Tr [aVV] = X> 


adj 


— 2n 0 . 


(7) 


( 8 ) 


(9) 


For A < 0 the ground state is doubly degenerate, i.e. no = ±1. (For the SU(3) model (Eqn. 
[6]) this degeneracy corresponds to two possible choices of the g matrix: g = e ± 2 W 3 /.) For a given 
choice of the sign one can consider deviations of the field g from the vacuum configuration as small. 
Then the low energy theory becomes a theory of three weakly interacting bosons with a Lagrangian 
density given by 


C eff = -^(<9 M n) 2 + 2|A|n 2 + 


( 10 ) 


where the dots stand for higher order terms. By rescaling n a —>• n a /k 1 ^ 2 we see that these terms 
contain powers of A:” 1 <C 1. Here the mass scale, 


m ~ 



|A[ 

k 


( 11 ) 


is obviously the one which is envisaged by RG considerations as k —>• oo (see Eqn. |4}. The double 
degeneracy of the vacuum is in agreement with the structure of the potential in the abelian action 
(Eqn. [5}, which is periodic under translations with a two-dimensional lattice generated by the 
vectors eW with two minima in each elementary cell, as shown in Fig. [U 

For A > 0 the situation becomes more interesting. At energies below m, the field uq is suppressed 


and n 
one 


recomes a unit vector. As a consequence the Wess-Zumino term ([2} becomes a topological 


28 |, 




T(i(j a n a ) = ^ fd’xe^n x = hr©. 


( 12 ) 


The coefficient of the topological term in Jl} is nk. Since the value of 0 is always integer, iA:7T0 
contributes nontrivially to the action only when k is odd. In summary, the low energy effective 
action for A > 0 is 


S = J d 2 x(d II n) 2 + brA;©, n 2 = 1. 


(13) 
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FIG. 1: The potential surface for k = 4. The black lines show the elementary cell of the periodic potential, 
blue dots show the two minima relevant for A < 0 that are connected via a saddle point. The red dot is a 
maximum of the potential, which becomes the only minimum (per elementary cell) for A > 0. 


This model is exactly solvable. For k even 


311 ] the particles are massive triplets with mass 


M tr ~ mfcexp(— k/2), m ~ A 1 / 2 . 


(14) 


The triplet structure of the particle multiplet agrees with the result for k = 2. However note that 
the mass scale at A > 0 is much smaller than the RG scale (Eqn. 0|). For k odd, the scale m (1141) 
marks instead a crossover into a basin of attraction of the critical point of the <S77(2)i WZNW 


model 


32]. 


III. TCSA FOR THE PERTURBED WZNW MODEL 
A. The truncated conformal space approach for SU( 2)*, perturbed by the trace of adjoint 


For a numerical determination of the spectrum we use the truncated conformal space approach 
, adapted to the SU( 2)*. WZNW model. For perturbations of WZNW models with levels k = 1 
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and 2 TCSA was applied previously in [33[; for the present work we developed a general purpose 
TCSA code working for all k and any perturbing operator. On a Euclidean space-time cylinder of 
circumference, R, in the spatial direction, x, the Hamiltonian has the form, 


fR 


H = Hk + A / dx&( 0, x), 


(15) 


where H *. is the Hamiltonian of the SU( 2)*. WZNW model (Eqn. [2]) and the perturbing operator 
<3? is minus the trace of the adjoint field (Eqn. [3]). Due to translation invariance, the Hamiltonian 
is block-diagonal on eigenspaces of the conformal spin Lq — Lq. The symmetry algebra is generated 
by Kac-Moody currents, J a (z), which satisfy the OPE, 

f a/3 


j^ {w ) = k ^. + fU2M + 0 ( i). 

2 (z — w) z z — w 


(16) 


where q a>3 is the invariant metric of the Lie-algebra su( 2), and are the structure constants. In 
the basis, the sit(2) algebra relations can be written as 


[S a , s?] = ffs\ f? = 0 a + p + r, 


(17) 


f°+ + = -/r = -/r = /r u = 1; /<T = -fo + = 2, 


c+0 


f 0- *-0 


and the metric is 


q a(i = 0 a + f. 3/0; 

q 00 = 1; q ±T = 2. 


(18) 


The modes of the current obey the Kac-Moody algebra, 

W = f £(c - *)”■/”( 0; 


(19) 




= ffjn+m + -"K^n+m.O, 


(where 2 denotes an arbitrary reference point of insertion; whenever omitted, it is taken to be 
2 = 0). The energy-momentum tensor is given by the Sugawara construction, 


1 


r 0) = j—^Qap ■ J a (z)J'\z ) 
d( 


( 20 ) 


L 


2(/ri 


where the modes L n satisfy the Virasoro algebra, 


k + 2 


jP 

'n—m •? 


raEZ 


— -A'n+m “1“ ^2^^ r; 


3fc 

/c + 2 ’ 


c = 


( 21 ) 
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We recall that for level k integer there are k — 1 primary field multiplets z), with j = 0, 

1/2, .. k/2 and m,fh = — j , —j + 1,..., j. They are normalized as 


( Zl > (*2, Z 2 )) 

Ki) 


^jlj2^mirri2 $m\m2 


(zi2Zi2)~ 2hUl) ; 


j{j +1) , 

k + 2 ’ 


( 22 ) 


The Hilbert space of the conformal field theory is given by 

®V,XV r (23) 

i 

where V ; - is the irreducible representation of su(2)k with highest weight j, and the bar denotes 
the antiholomorphic component. The ground state level of the module Vj (g> Vj is spanned by the 
multiplet, 

I m,fh)j = <h^(0,0)|0), (24) 


where |0) is the SL(2,C) invariant conformal vacuum state. The module is generated by the 
raising operators m < 0; care must be taken to factor out null vectors to obtain the irreducible 
representation. 

We remark that for computational simplicity the left and right Kac-Moody algebras were im¬ 
plemented in the same way. This is in contrast with the usual WZNW formalism, where the 
fundamental field transforms as 


9 9l99r ■ 


(25) 


However we take our fields to transform as 




(i) 


D {j) ( g L )mm ' DW (g R )rnfn'$ { £,rn '> 


(26) 


where D^ 1 is the SU(2) representation corresponding to spin j. This is related to the usual 
definition by a contragredient transformation applied to g R which is equivalent to a redefinition of 
the basis for SU(2). 

n 

The operator product coefficients of primary fields were derived in }34j . Introducing the gener¬ 
ating function fields 


*W(x,x;z,z) = '£ t J( 2 [ ) ( 2J \xi-™xi-^%{z,z\ 
“ V \ m + JJ ym + J/ 

m.m V x ° 7 v 17 7 


( 27 ) 
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which have the two-point function, 

^ n \x 1 ,X 1 -,Z 1 ,Z 1 )^ :l2 \x2,X2-,Z2,Z2)^ = SP (x 12 X 12 f n {z 12 Z \ 2 )~ 2h{:,l) , (28) 

the operator algebra is fully specified by the three-point functions of the primary fields given by 

^ n \xi,xi, zi, zi)^ n \x 2 ,x 2 , z 2 ,z 2 )^ j3 \x 3 , x 3 ,z 3 , z 3 
= C(ji, j 2 , j 3 ){x 12 xi 2 ) jl+j2 ~ j3 (xi 3 xi 3 ) jl+j3 ~ j2 (x 23 x 23 ) j2+j3 ~ jl 

x(^122l2) /l03) " /l0l) " /l02) (-2l3^13) /l02) “ /l(il) “ /l03) (^23^3)' l(jl)_Mj2) ' ,l03) , (29) 


where the structure constants C are given by 

1 


C(j 1, j'2 ; h) = 7 


k + 2 


T}/ ■ I • I • I u 2TT ^ > 0' 1 +i 2 +i 3_ 2jn) 2 

P{n+j 2 +js + i) J_1— , lN -; 

«»1 7(%i5 1 )^(2i„) 2 


l(x) = 


r(x) 


r (i - x) 


p{j)= n 7 


n= 1 


n 


fc + 2 


(30) 


and are fully symmetric in their arguments. Structure constants for component fields can be 
obtained by expanding in the variables x, x and using Eqn. [271 

The trace of the adjoint field can be expressed as the component of the j = 1 field which is a 
singlet under the global SU(2) symmetry generated by Jq + Jq : 


1 


<T) — _22_ ( (f)! 1 ) _l_ (T)^ _ <T>^ ) 

^ + “hi ^0,0 J i 


(1) 


s(i) 


(31) 


where the prefactor ensures that the conformal two-point function of the perturbing field is canon¬ 
ically normalized: 


(<h(^ z)*(w, w)) = h = . (32) 

In fact, the field 3> as defined in Eqn. [31] differs from the trace adjoint defined previously by a sign, 
which is consistent with the form of the Hamiltonian in Eqn. [15j This sign can be identified from 
matching the TCSA result against the semiclassical predictions, performed in the sequel. 

Using complex coordinates £ = r + ix, after the exponential mapping from the plane to the 
cylinder, 


* = e 2 </ R , 


(33) 


the Hamiltonian can be written as 


H = 


2vr 
"R 


Lq + Lo — 


27r R 


2-2 h 


4(fc + 2)J + X R (2 7 r) 1 - 2fe$(1,1) ' 


(34) 
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The dimensionful coupling constant A can be used to define a mass scale M 

|A| = M 2 ~ 2h . (35) 

In all our subsequent computations we use dimensionless quantities measured in units of M. The 
dimensionless volume parameter is given by r = MR and the Hamiltonian can be written as 


H 

M 


2vr 
r 


„2-2 h 


L ° + L ° ~ 4(fc + 2) + Slgn(A) (2 t r) 1 -^ $(1 ’ X) J ' (36) 

The matrix elements of the perturbing operator between descendant states can be computed by a 
recursive procedure using the relations given by the Kac-Moody algebra. Truncating the Hilbert 
space at some descendant level N, the dimensionless Hamiltonian becomes a finite numerical matrix 
for any given value of r, which can be diagonalized numerically, resulting in a raw TCSA spectrum 
that depends on the truncation level. In many cases the raw TCSA data already gives an accurate 
spectrum; however, we used both numerical and analytic schemes to eliminate cut-off dependence 
and obtain better results. 

Since the perturbation is a singlet, it conserves the z component of the diagonal SU(2) (gi = gn) 
which is 


Q = J§ + j§. (37) 

All eigenvalues of the charge Q are integers and the Hilbert space can be decomposed into sectors 
labeled by the eigenvalues of Q. In all our computations we only show results for states with Q = 0. 
Since the Hilbert space decomposes into integer spin representations of the diagonal SU(2), each 
such multiplet has a single level lying in the Q = 0 subspace. It is guaranteed analytically, but we 
also checked numerically that any Q / 0 state is degenerate with one of the levels in the Q = 0 
subspace, and that degenerate states form full multiplets. In addition, all TCSA data we present 
correspond to zero-momentum states (i.e. Lq — Lq = 0), as non-zero momentum subspaces contain 
no new physics. 

It is also important to observe that according to the Kac-Moody fusion rules under the per¬ 
turbation (Eqn. [3ID , the Hilbert space (Eqn. l23l) decomposes into even and odd sectors which 
originate from states with j integer and half-integer, respectively. This will be important in the 
sequel. 


B. Testing the TCSA 

In this subsection we describe a number of tests that we put our TCSA code through. Such tests 
are very important as there are no exact results to verify the TCSA due to the non-integrability 
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of the theory. The first test that we considered is specific to k = 4. As we have discussed, we have 
two realizations of SU(2)±+T?<& a( ij: one treating the conformal basis of SU( 2)4 in the language of 
current algebras (the picture we are focusing on in this paper) and one treating this same basis as 
a two boson theory, one boson with compactification radius of \Z2ir and one orbifolded boson with 
compactification radius of 2\/6it (see Eqn. E}. While these are very different starting points for 
the TCSA, they must lead to the same answer. And we have checked, at least in the sector of the 
theory containing the ground state, that they do. While this is an important check of the code, it 
only applies to k = 4. We have thus also considered the case SU(2)i perturbed with the singlet 
component of the j = 1/2 primary field (equivalent to sine-Gordon theory in the SU(2) symmetric 
point of the attractive regime) as well as the k = 2 case, where the spectrum must agree with that 
of three Majorana fermions; some details on this latter case are given in subsection IIV Dl Both of 
these tests indicate the code is working. Finally, as discussed in subsection IIII C 41 we show that 
corrections to the ground state due to changing the cutoff in the theory as determined by the TCSA 
match those computed analytically using conformal field theory. This analytical computation is 
non-trivial and so provides another important check on whether the code is behaving as it should. 

C. Renormalization methods 

The first improvement to the raw TCSA is given by the numerical renormalization group (NRG) 
method introduced in jf|. The procedure consists of starting at a cut-off value where the whole 
matrix can be diagonalized and then incorporating higher energy levels in chunks of a given step 
size (number of states added at each step) until the target value of the cut-off is reached. This is 
necessary as the number of states grows very fast with the cut-off. For example, in the integer j 
sector of the k = 4 theory for descendant levels N = 3,4, 5, 6 , 7 we have 1427, 6373, 23498, 83144 
and 264129 zero-momentum states respectively in the Q = 0 sector. Our computational capacity 
allowed us to reach the descendant level N = 5 with exact diagonalization, while to reach N = 6,7 
we have used the NRG procedure. 

The next improvement takes into account the contribution of states above the cut-off using 
perturbation theory to second order in A. There are several schemes in the literature that can be 
used, but a compromise must be struck between computational costs and accuracy. Below we give 
a short description of each procedure, and compare them in order to make an optimal choice for 
our problem. 
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1. Vacuum counter term 


The contribution from the omitted high-energy states dependence of the ground state energy 
was computed using the results in [yj], and the counter term necessary to eliminate the cut-off 
dependence to second order is given by 
irR 3 ~ 4h 


5E 0 = 


irR 3 ~ 4h 1 (V{2h + N + 1) V 

2(27r) 2 ~ 4ft h + N + 1 \T{2h)T(N + 2)) 


x 4 F 3 ( 1,1 + h + N,l + 2h + N,l + 2h + N;2 + N,2 + N,2 + h + N;l) 
{2n) 4h - l R 3 - ih 


-N' 


Ah—2 


+ ... 


4(2 h - l)T(2/i) 2 
with 4 T 3 denoting a generalized hypergeometric function. 


(38) 


2. Counter terms for excited states 

To eliminate cut-off dependence for excited states we can use a scheme developed in [15|]. The 
idea is to separate the Hilbert space into a low-energy part (labeled by l ), which is included in 
TCSA, and a high-energy part (labeled by h) which consists of states above the cut-off. For any 
state we split its eigenvector c into low- and high-energy parts q and c^; similarly, the Hamiltonian 
can be split into a block form according to 

Hu Hi h 

Hhi H hh 

The full eigenvalue problem can be split accordingly 


H = 


(39) 


Huci + Hi h Ch = eci, H h ici + H hh Ch = £Ch, (40) 

where e is the exact eigenvalue. Eliminating the high-energy components, c^, gives 

[Hu —H lh (H hh - ey 1 H h j]c, = eq. (41) 


We write 


A-Hfuii 


H = H 0 + V, 


(42) 


where Hq is the conformal Hamiltonian and V is the matrix of the perturbing field. Note that the 
off-diagonal components only involve the perturbing matrix, so we obtain 

AHfuii = —Vi h (H 0 + Vh h — e) _1 14; 

~ —Vi h (Ho — s)- 1 V hl +0(E^ X ), 

" -V-' 

AH 


( 43 ) 
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where E max is the cut-off in energy units, and we used the fact that V b h only contributes at higher 
order. First order perturbation theory in E^ ax gives 


£ = EtCSA + ctcsa A H ctcsa > 


(44) 


where 


HilCTCSA = EtcSACtCSA- 


(45) 


To calculate AH we can approximate e -A Etcsa in Eqn. 03] 

M(E) ab 


A H ab = - 


j 

JE, 


E - Etcsa ’ 


where M{E) ab is defined by 


POO 

(a\V(r)V(0)\b) = / dEe-( E ~ E ^ T M(E) ab . 

J o 


'0 

Here V (r) is understood as a Euclidean time-evolved version of the perturbing operator, 

r R 


Then one can write 


H(0) = f dx$(ix). 

Jo 

pR pR 

V(t)V( 0) = / dx i / dxoQ^ixi + r)<&(ix2), 
Jo Jo 


(46) 


(47) 


(48) 


(49) 


where 4* has scaling dimensions, [h, h ). The leading contributions can be computed by considering 
the most singular terms in the operator product expansion on the cylinder 


4>(ixi + r)4>(ia: 2 ) ~ \ix\ - 1 x 2 + r | ih+2hv ip(ix 2 ) + ... (50) 

where ip has scaling dimensions and we only consider the first few operators. Indeed, the 

leading contribution is the identity with h\ = 0 and 6 * 1 $$ = 1 ; since we neglect the subleading 
terms of the identity contribution, we cannot include any operators with h v > 1 / 2 . 

In the spin 0 sector, using translation invariance of the states |a) and | b) gives 


(a\y{ix 2 )\b) = (a|</>( 0 )| 6 ) 


(51) 


2tt\ 


2 hu 


— ( r J ( a | ^Pplane ( 1 ) 1 ^)- 


The integrals can be explicitly computed 
r R r R 


r„ xi r 

Jo Jo 


dx 2 \ix\ — ix 2 + r|“ = 2 R 2 r a 2 F 1 ( -, — —; - 5 - ) + 


1 a 3 R 2 


2’ 2’ 2’ t 2 ) ' a + 1 


r 2cx.-\-2 


(f ? 2 + r 2 ) a+1 


OL -hi 


r f— 

= \pK R - 7 2 / r l+a _|_ J egg gi n g u l ar _ 

r(-f) 


( 52 ) 
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To leading order we can neglect the term e EaT in the correlator, since the energies E we consider 
are above the cut-off E max , while E a is below it; similarly we can take Etcsa —> 0 in Eqn. 1461 
This gives 


A ^ _ (2ir) 4h - 1 R 3 ~ Ah N 4h - 2h *- 2 

ab ~ 4(2 h - 2/v - l)T(2/i - 2h v ) 2 


ane (m, 


(53) 


for the contribution coming from the operator tp in the OPE. Note that this equation gives the 
counter term in an operator form, in an arbitrary basis of states |a). The eventual correction to 
any given energy level is computed by evaluating its matrix in the TCSA basis and computing its 
expectation value with the TCSA eigenvectors corresponding to the given level at cutoff N. Note 
that when applied to the case of the identity, this result reproduces the leading behaviour of the 
vacuum counter term (Eqn. 


3. Running coupling 


To the leading order, the counter terms (Eqn. l53l) are independent of the state under consider¬ 
ation and can be considered as local operators added to the Hamiltonian. This forms the basis of a 
renormalization group (RG) by defining a level dependent running coupling such that the counter 
term describing the cut-off dependence is compensated by changing the coupling. Let us denote 


h = 


t) 


2h—2 


(54) 


and consider the case when (/? = $, i.e. = h. Then we can introduce a cut-off dependent 
coupling by requiring that the contribution of high-energy states from level N is compensated by 
the change in the coupling. The contribution from the iVth level to the <h term in the Hamiltonian 
is just the derivative of Eqn. [53] with respect to N (to leading order in 1 /N), which is equal to 


A/fJ = (1)14). 


(55) 


We make the coupling depend on N and stipulate that it is evolved from N to N — 1 by including 


the counter term’s contribution. This leads to the RG equation of the form used 

, ,.2 i\j2h—3 , 

Rn — Rn-i + Hn -i Y(h ) 2 ^ +••• 

where the dots denote terms subleading for large N, which can be integrated to 

Hn ‘2'kCqqq 


in [E, 12, 13] 


Moo — 


1 + Ci^vA^" 2 ’ 


Ci = 


{Ah — 4)T(/i) 2 ’ 


(56) 


(57) 
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to leading order for large N. Let us denote the dimensionless energy levels by ej(r), which are 
just the eigenvalues of the dimensionless TCSA Hamiltonian (Eqn. l36l) as functions of r, with the 
vacuum being eo (r). 

This RG equation, inasmuch as it can be derived from the invariance of the partition function 
under changes in coupling 


12j, expresses the RG invariance of the gaps of the form: 

EiiXoo) =eI N) (X n ). 


(58) 


Due to Eqn. [54l this invariance can be reinterpreted in terms of the dimensionless energy levels 
via 


ei(roo) = —e\ N \r N ), 

T oo 


(59) 


where e* are the energy levels at cut-off N = oo. Note that the energy also needs to be rescaled 
due to the 1/r prefactor in the dimensionless TCSA Hamiltonian (Eqn. 1361) . 

We remark that at higher orders in 1/N the counter terms (and therefore the running coupling 
as well) are state dependent [h|. One way to take this into account is to compute the full counter 
term (Eqn. ITTIi without using the approximation of the previous subsection, i.e. keeping the 
dependence on Etcsa and E a in Eqns. 06] and 07] This leads to a rather complicated and 
computationally expensive method even for the counter terms themselves, and it makes rather 
difficult the implementation and solution of the corresponding renormalization group equations, 
which describe non-local Hamiltonian terms [ljj. Henceforth we neglect these higher corrections 
in our computations. 


4- Renormalizing the ground state 

The ground state of the theory has the conformal vacuum as its ultraviolet limit, and is contained 
in the even, Q = 0, zero-momentum sector. In Fig. 0]we show results coming from NRG+TCSA, 
the effect of the vacuum counter term (Eqn. 1381) . as well as the results obtained by implementing 
both the counter term and the running coupling according to Eqn. [59] In small volume we can see 
that the counter terms at different cut-off levels scale the energy level to the same curve, verifying 
that the subtraction provides reliable results even when starting from NRG+TCSA data with low 
values of the cut-off. Note that taking into account the running coupling gives a further significant 
reduction of cut-off dependence. This scaling is one of the verification tools we used to affirm our 
belief that our code is producing correct results. 
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FIG. 2: k = 3 (top) and 4 (bottom), A > 0 ground state energy data obtained by TCSA+NRG for cutoff 
levels N = 5,6, 7,8 and N = 4, 5, 6, 7 (dashed lines); the results with the counter term (Eqn. [38]) (red lines); 
and the data involving the counter term and the RG improvement (Eqn. 1591) (green lines). The insets show 
the same results blown up on the 1 < L < 14 interval. Note that in the direction of increasing cut-off N the 
subtracted (red), and the subtracted and renormalized (green) levels move less as the cut-off grows, which 
is a further confirmation of the validity of the renormalized TCSA. 

E/M 




So 

A > 0 

A < 0 

k = 3 

-4.55 ± 0.01 

-7.69 ± 0.01 

k = 4 

-2.21 ± 0.01 

-5.03 ± 0.01 


TABLE I: The bulk energy density £q given in units of M (cf. Eqn. l35l) . 


The slope extracted from the linear regime of the vacuum level gives the ground state (bulk) 
energy density £q, and can be estimated by fitting a linear function in the appropriate range of 
volume. The resulting estimates are given in Table [H 
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5. Numerical application of the renormalization methods 

To eliminate the additive bulk energy renormalization (Eqn. 1381) we consider the gaps relative 
to the vacuum. The running coupling (Eqn. l57l) leads to the renormalization prescription 

ej(roo) - e 0 (roo) = — (e[ N \r N ) ~ e^\r N ) \ , (60) 

r oo ' ' 

which we call the RC (running coupling) correction. It is also possible instead to add the counter 
terms (Eqn. l53l) where we can include the contribution of all operators ip below a certain h^ 
chosen to keep the slowest-decaying contributions (in practice we chose to incorporate the primary 
contributions). This will be called the CT (counter term) correction. The difference between the 
two corrections is that in contrast to the CT correction, the RC correction only involves a single 
operator contribution, however by introducing the running coupling it sums up the leading power 
in the cut-off dependence to any order. 

We illustrate the renormalization method for the first excited level for k = 4 and A > 0, which 
comes from the odd sector and consists of three degenerate states with Q = +1,0,—1 forming a 
triplet under diagonal SU(2). We consider the state corresponding to a stationary particle, which 
can be found in the zero-momentum sector. For higher cut-offs we find that the two prescriptions 
converge to each other as illustrated in Fig. [3l so we can use the computationally simpler RC 
method to obtain renormalized results. 

IV. NUMERICAL RESULTS FOR THE SPECTRUM 

A. The A < 0 case 

In this regime we expect two triplets of weakly interacting particles that have a mass given by 
Eqn. [TTJ We show the spectra of zero-momentum Q = 0 states relative to the absolute ground 
state in Fig. [4] which clearly show a doubly degenerate vacuum structure. In finite volume, the 
degeneracy of the vacua is lifted by the tunneling, which vanishes exponentially with the volume. 
Due to the Z 2 symmetry relating the two vacua |1) and 12), the finite volume ground states are 
given by 

|±) = + (|1>±|2», (61) 

and are expected to emerge from the even and odd sectors, respectively. Since we plot the energies 
relative to the absolute finite volume ground state |+), the presence of these vacua is signalled by 
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FIG. 3: Finite volume gap in the perturbed SU{2)^ model for A > 0 at cut-offs IV = 2, 3, 4 and 5. The plot 
shows the raw TCSA data (without NRG, blue circles), and those after the RC (dashed lines) and the CT 
(magenta squares) corrections were applied. The gap can be estimated to be A = 0.36 ± 0.03, and the error 
is approximated by the difference between the gap estimates for N = 3 and 6. 


a state originating from the odd sector with a relative energy approaching zero exponentially with 
the volume. 

As expected from the semiclassical considerations, the first excitations are indeed a triplet of 
particles, and they too appear in two copies according to the vacua. Their triplet nature can be 
seen both by looking in the Q = ±1 sectors for the other components of the multiplets, but also 
from the fact that the energy levels in the ultraviolet (mR ~ 0) limit are seen to emerge from 
conformal states transforming as a triplet under the diagonal SU(2). In particular, the lowest 
lying states come from a quartet of states created by the primary field and a nonet of states 

created by the primary field <J>^. Under the global SU(2), the quartet decomposes as 

I®I = 0©1, (62) 

with the singlet giving the second vacuum state, while the triplet corresponds to the first triplet of 
one-particle states. In the plots of Fig. []]the quartet corresponds to the first two blue lines (their 
colour indicates that they come from a sector created by a primary Held with half-integer j). 
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FIG. 4: Finite volume spectra in the perturbed SU (2)*, models with k = 3,4, 5 for negative coupling constant 
at cut-offs n m ax = 7, 6, 5, respectively. We show raw TCSA data since renormalization here had an effect 
that is not visible on these figures. Colors represent energy levels in the integer (black) and half integer (blue) 
sectors. The red arrows show the gaps corresponding to one-particle states. For k = 5, two of these are 
already higher than the two-particle threshold (shown as the thick dashed line), and due to non-integrability 
they are expected to correspond to resonances. 
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A 



FIG. 5: TCSA mass gap for negative coupling constant as a function of fc _1//2 . We show data coming from 
k = 2, 3, 4 and 5 models. We also put error bars on the data points which we calculated by subtracting 
the gap estimates with and without applying the RC improvement, but they are so small that they are 
practically invisible in the plot. 

The nonet 


1 0 1 = 0 © 1 © 2, (63) 

contains the other triplet of one-particle levels while the single and quintet are excited states. 
They correspond to the lowest three black levels visible in the plots, which are exactly (the Q = 0 
components of) the three multiplets. Of these three levels, the one corresponding to the triplet 1 
in Eqn. [63] approaches the other (blue) triplet from Eqn. 1621 exponentially with increasing volume. 
They also level out exponentially, signaling that these are single-particle states, coming in two 
copies according to the degenerate vacua. 

In Fig. [5] we show that for larger values of k the gap measured from the flat portion of the 
first one-particle levels indeed follows the fc -1 / 2 scaling of the particle mass expected from the 
semiclassical arguments. The spectra also show the presence of additional states below the two- 
particle threshold. Due to the double degenerate vacua, one expects kinks interpolating between 
them. With the periodic boundary conditions imposed by TCSA one can only see states with an 
even number of kinks such that the sequence of vacua interpolated by them has the same starting 
and end points. In addition, these kinks are also expected to have bound states, and the lowest 
lying particle triplet can be identified with the lowest mass kink-antikink bound states. 

In the absence of more detailed knowledge about the theory, at present we cannot identify the 
higher states conclusively, but it seems that at least the first few levels very much resemble the 
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breather doublets seen in the 2-folded sine-Gordon theory 


351 ], so it is likely that these are indeed 


higher kink-antikink bound states beyond the lowest triplet. These levels can be seen in Fig. [4] 
as pairs of black and blue lines approaching each other and also leveling out exponentially at the 
same time. The fact that one of these always comes from the even, while the other from the odd 
sector (as shown by their colours) confirms the interpretation that these are indeed two copies of 
higher kink-antikink bound state multiplets (note that these are not necessary triplets; from the 
identification of the nonet lines we know the next two are a singlet and a quintuplet). 

One can also see that the number of such one-particle level candidates increases with k, which 
is what is expected in the semiclassical limit 3(], 37]. In addition, the characteristic dependence of 
the lowest particle mass on k suggests that the mass scale M is related to the kink mass, and that 
the spectrum of bound states becomes dense for large k, analogously to the <h 4 and sine-Gordon 
models treated in [3a, [37|]. Due to non-integrability of the model it is also expected that two-kink 


bound states over the two-particle threshold are in fact resonances whose finite volume signatures 
must resemble those in the two-frequency sine-Gordon theory studied in 38|; however, our data 


do not allow a reliable identification of these signatures at present. We also remark that in spite of 
non-integrability, there are also apparent level crossings in the spectra, e.g. between even (black) 
and odd (blue) levels since they do not mix under the perturbation. The additional higher states 
can be interpreted as two- and more particle levels composed of particles and/or even number of 
kinks. 


B. A > 0, even k 

In Fig. [6] we show the results for positive A and even level k. We observe a single vacuum 
and a triplet of one-particle levels, which are consistent with the effective cr-model picture based 
on the action in Eqn. [13j From the first excited levels alone, the gaps can be estimated as 
A = 0.37 ± 0.03 and A = 0.18 ± 0.02 for k = 4 and 6, respectively, which are much smaller than 
those for A < 0 and decrease strongly with increasing k as expected from the semiclassical result 
(Eqn. fl4l) . Unfortunately, the numerical accuracy for higher levels is not very good, which at this 
stage precludes their interpretation; this is further complicated by the smallness of the gaps, which 
means that the volumes we could reach are in fact very small in terms of the correlation length, 
therefore sizable exponential finite size effects are expected. 

For k = 4, the two-particle gap can be seen to be of order 0.9 from the lowest level above the 
one-particle state; this is roughly consistent with the gap estimate above, but cannot be trusted 
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FIG. 6: Spectra in the k = 4 and 6 models for A > 0 at cut-offs N = 7 and 4. Data was improved by the 
RC method. Colors represent level data obtained in the integer (black) and half integer (blue) sectors. 

to be of the same precision as there are clear signals of large residual cut-off dependence in the 
behaviour of the level, e.g. the fact that it curves upwards for larger values of mR, while in reality a 
two-particle level must approach the threshold from above with a behaviour ( mR )~ 2 . For the case 
k = 6, the truncation achieved is very low and so the higher levels cannot be trusted, precluding 
their analysis for the time being. 


C. A > 0, odd k 

In this regime, as we discussed in Section II, we expect a massless flow to an SU( 2)i low-energy 
fixed point. The mass scale (Eqn. [DO) in this case corresponds to the cross-over scale. The spectra 
shown in Fig. [7] do show marked differences from the even k case. All levels are monotonically 
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decreasing with the volume, rather than leveling out as in a massive spectrum. Also, one would 
expect the gap for k = 3, 5 larger than for k = 4, 6, respectively, while compared to the data in 
Fig. [6] one readily sees that the distance between the ground state and the first excited state is, 
instead, markedly smaller and monotonically decreasing even at the largest volume shown. 

The existence of an infrared fixed point implies that for large values of the volume the energy 
levels should behave as 


Ei — En 


2irxi 

~R~ 


+ 


where Xi are scaling dimensions in the SU{ 2)i theory, and the dots indicate corrections to the 
low-energy scaling limit. One can define the scaling functions 

Ft 

Di = —{Ei -E 0 ) = Xi + ... 

Z7T 

As shown in Fig. [HJ the detailed matching is rather limited. There are reasons for which this is 
expected. First, from the gaps measured for even k, the typical scale parameter for the cross-over is 
expected to be MR > 4. In addition, the cross-over itself is slow, due to the fact that the irrelevant 
perturbation describing the incoming direction in the infrared is the current-current perturbation 
of SU( 2)i, which is only marginally irrelevant and leads to a logarithmic approach to the fixed 


point 


401 ]. Therefore one expects that the fixed point would only be observable for volume values 


much higher than allowed by TCSA accuracy. 

It has long been known that observing an infrared fixed point in TCSA is difficult 
Ref. [39f] an attempt was made to observe the flow from the tricritical Ising conformal minimal 


39}]. In 


model to the Ising minimal model by perturbing the tricritical Ising theory with the subleading 
energy perturbation, with the conclusion that the behaviour of the first excited state was not 
inconsistent with the existence of the fixed point. In a study of two-frequency sine-Gordon model 
4j], the Ising fixed point was just barely in the reach of the TCSA. In that case however, the 
sine-Gordon frequency provided a parameter which could be tweaked to improve convergence to 
the point that the first two scaling dimensions of the infrared fixed point could be extracted, albeit 
with considerable errors. Note also that in these examples the approach to the fixed point was 
much faster (given by power corrections). 

Pending more accurate TCSA numerics (which would require a more accurate modeling of the 
cut-off dependence, and more computing power to allow higher truncation levels), we can only 
say that the TCSA data are qualitatively consistent with the existence of a low-energy quantum 
critical point, and the first scaling function in Fig. [8] is also roughly consistent with the lowest 
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FIG. 7: Spectra for k = 3, 5 for positive coupling constant at cut-offs N = 8 and 6. Data was improved by 
the RC method. Colors represent level data obtained in the integer (black) and half integer (blue) sectors. 


scaling weight 



D. The case k = 2 

For this case one expects a spectrum of free massive Majorana fermions; in our units given by 
Eqn. [35] the fermions mass is 

27T 

m = —=M. (64) 

V3 

The resulting spectra for A > 0 and A < 0 are shown in Fig. [9j They are exactly the spectra 
expected for three Majorana fermions in the Z 2 symmetric and Z 2 symmetry breaking phases, 
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D 




FIG. 8: Same data as in Fig. [7] but here we show the scaling functions, D, = (E t — Eq). The dashed 
lines are the first few scaling weights in the SU{ 2)i CFT. 

respectively. Note that in the A < 0 the excitations are kinks, therefore there are no single-particle 
levels and every multi-kink level appears in two copies (one even, the other odd), which are split 
by tunneling effects decaying exponentially with volume. There are no kink-antikink bound states 
below the two-particle threshold. For the A < 0, the k = 2 point is analogous to the free-fermion 
point of sine-Gordon theory, while the k > 2 cases correspond to attractive regime. The main 
difference from the sine-Gordon case is the non-integrability and the presence of SU{2) invariance. 

In the symmetric phase A > 0, there is a unique vacuum and the even/odd sectors correspond 
to levels with even/odd fermion numbers. In particular, the odd sector does contain one-particle 


states. 
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FIG. 9: Spectra for k = 2. The red dotted line (shown only for A > 0) corresponds to the fermion mass 
to = 2nM/y/3, while the red dashed lines show the position of the two-particle threshold 2 to. 


V. SUMMARY AND CONCLUSIONS 


We have studied the SU(2)k+Tr& a dj theory using the TCSA+RG approach. We have compared 
our numerical results to the semiclassical analysis of this model. We recall that the semiclassical 
picture suggests two regimes. One of them corresponds to the negative sign of the coupling and has 
a doubly degenerate ground state characterized by a nonzero vacuum average of the SU(2) matrix 
field (Tr g) = ±<r. Therefore one expects kinks interpolating between the different vacua. Indeed 
TCSA data show a spectrum that resembles those in the <h 4 and also the 2-folded sine-Gordon 
theory, where particles arise as bound states of kinks with a spectrum that becomes dense in the 
semiclassical limit which corresponds here to large k. The lowest lying excited states consist of 
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two triplets of particles with the same mass and can be identified with the lightest kink-antikink 
bound states. 

In the other regime for even k the semiclassical considerations suggest a triplet of massive 
excitations, whose low-energy dynamics is governed by the 0(3) sigma model. The presence of 
a single vacuum with a massive triplet is confirmed by TCSA, but the precision is significantly 
smaller since due to the much smaller gap, it is necessary to go to much larger volumes which 
increases the truncation effects. Even so, the vacuum energy density and the mass gap can still be 
extracted with a reasonable accuracy by applying renormalization group improvement techniques. 
For odd k, the semiclassical considerations imply the presence of a low-energy quantum critical 
point described by the SU(2)i conformal field theory. The TCSA data are consistent with this 
prediction, but are not accurate enough to identify the nature of the fixed point conclusively. We 
still think that the semiclassical picture has proven robust enough so that this prediction can be 
trusted. 

For the exactly solvable case k = 2 the spectrum of the model describes three Majorana fermions 
which also constitute a triplet. The model is then non-interacting, therefore there are no more 
particles in the spectrum. For the two signs of the coupling the spectrum differs as usual for a free 
fermion theory, and fits well into the pattern observed for k > 2. 


Appendix A: Realization of SU(2)k + Tr$ a dj 


Here we give an example of how the perturbed WZNW model of Eqn. [T] can appear in the 
context of an electronic model. We start with a lattice model with U(l)xSU(2)xSU(k) symmetry, 


H = 


J2{~ t + !)Vv( n ) + H - c ■ 

n 

- J ip) a {n)ip ia (: n)il)\ a , (n)ip ja > (n) j, 


+ U [^]a{ n )^A n )][^\a'{ n )^'in)] 

{jo- }^{ier'} 


(Al) 


where tpj a (n) and are creation and annihilation operators of electrons located at sites n; 

a = ±1/2 are spin and i,j = 1 ,...k are orbital indices. Treating the interaction as small in 
comparison with the Fermi energy and assuming that the band is far from being half filled, we 
separate fast and slow Fourier harmonics of the electron operators: 


ip(n) = e~ ikFnao R(x ) + e ikFnao L(x), x = na 0 , 


(A2) 


where kp is the Fermi wave vector and a q is the lattice constant, and arrive to the continuum version 
of Eqn. IA1I in the form of the chiral Gross-Neveu model with the most general current-current 
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interaction. The corresponding Hamiltonian density is 


T~L — iv(—RpjdxRaj + d x L a j ) + + (A3) 

g 0 R t(r a 8) T) J R J L t (r a <g> J)L + 5so i?t( T a <g> a b )RL f (r a <g> a b )L + g s R t(/ <g> (/ ® a a )L, 

where cr a (a = 1,2,3) acting on the Greek indices and r a (a = 1, ...k 2 — 1) acting on the Latin 
indices are generators of the su(2) and su(k) algebras respectively, normalized as 

Tr(a a a b ) = Tr{T a r b ) = 


The coupling constants 31 , 2,3 are related to U and J while v the Fermi velocity is given by v = 
2 t sin(&Fao). 

Model (Eqn. IA3D is integrable for g a = g so /2, g s = g so /k where in the case the symmetry expands 
to U(l)xSU(2k). In this case the abelian sector is massless and the non-abelian sector is massive 
if at least one of g Q or g s is positive and g so 7 ^ 0 and is massless otherwise. It is also integrable if 
g so = 0. For this last case the Hamiltonian density can be written as a sum of three independent 
Wess-Zumino-Novikov-Witten (WZNW) models perturbed by the current-current interactions: 


U = ■ 

r 2 tt 
.k + 2 


2vr 

kT2 


(:«: + : J a L J a L :) + g s J a R J a L 


+ )+9onn 


+ 


'77 

-k 


■ 3R3R '■ + : 3L3L : 


+ 9c3R3L , 


+ 


(A4) 

(A5) 

(A6) 


where Jr/ l , Fr/li an< 4 Jr/l are SU( 2)^, SUjkh, and 17(1) left/right Kac-Moody currents. Each 


perturbed WZNW model is exactly solvable 
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25]. 


We consider the case g s < 0 and small g so when the term mixing the spin and the orbital sectors 
in Eqn. IA3I can be considered as a perturbation around the SU{2)k WZNW critical point. As we 
shall demonstrate, this perturbation is always relevant and is given by the SU(2)k adjoint operator. 


The standard analytic approach to the models of type 


On this basis certain robust predictions have been made 
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jqn. 


starts with RG equations. 


271 ]. In particular it has been argued 


that at the lowest energies the largest possible symmetry is restored (in our case it would be 
U(l) xSU(2k)). We do note that the reliability of such approach hinges in part on weak coupling 
as the Gell-Mann-Low function is universal only at first loop. The first loop RG equations for the 
model in Eqn. 


are 


9o = kgl+3kg 2 S0 /A, g so = (k - 2)g 2 so + g so (kg s + 2g 0 ), g s = 2g 2 + 2(k - l/k)g 2 a . 


(A7) 
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The case we are interested in is g o (0) > 0,g s (0) < 0. Then at g so = 0 the current-current 
interaction in the SU(2) invariant sector scales to zero and this sector is gapless. On the other 
hand, the interaction in the SU(k) sector (the orbital one) scales to strong coupling and the 
excitations in this sector become massive. This occurs at the RG scale £ 0 ~ l/g o (0)k. At finite g so 
the corresponding term acts as a relevant perturbation. Assuming that g so remains the smallest 
in flowing to the scale governed by £ 0 , we extract from Eqn. IA7I its value at this scale to be: 


9so(^o) 


g so { 0 ) 

g o (0) + 2\g s (0)\/k' 


(AS) 


We will assume that |g so (£o)l ^ 1 and consider the spin-orbit current-current interaction as a 
perturbation. As it was demonstrated in [T?j], this perturbing operator is the trace of the primary 
field of the SU(2)k WZNW model where is the field belonging to the adjoint representa¬ 
tion. This argument is based on the observation that the marginally relevant term with g os having 
scaling dimension 2 can be represented as a product of conformal blocks of the SU(k )2 and SU(2)k 
primary fields in the adjoint representation. This suggestion is based on their scaling dimensions: 


2k 4 

^■[5^)2] = —, d adj [SU(2) k ] = —, (A9) 

In the vacuum of the perturbed SU{k )2 WZNW theory, only the Tr ( h" rf? has a nonzero average. 
This leads one to the conclusion that after the high energy degrees of freedom of this theory are 
integrated out, the local operator Tr<f> arfj [SU (2)*.] will emerge from the product of the corresponding 
conformal blocks. 

We thus estimate the coupling of the perturbation to be equal to 


g,o{So)( r n* aiii [SU(k)]). 


(A10) 


As is shown in the main text, the physics of the model in Eqn. Q] depends crucially on the sign 
of A. If we consider the WZNW action in Eqn. [T| as a descendant of the fermionic model (Eqn. 
EH), the sign is determined by the product of signs of g so and the vacuum average of the adjoint 
operator in the SU(k) sector (Eqn. IA10I) . The ground state of the SU{k )2 model perturbed by the 
current-current interaction is degenerate and hence the magnitude and the sign of A depend on the 
vacuum. This degeneracy is lifted by the interaction with the SU (2) sector. As a result the sign 
of the interaction (Eqn. IA10D must be chosen so as to minimize the ground state energy. As we 
have demonstrated (cf. Table Hj), for a given k the ground state energy is always lower for A < 0 
where the masses in the SU(2 ) sector are the largest and the lowest energy excitations consist of 
a massive triplet of particles, which appear in two copies due to the degenerate pair of vacua. 
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